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Abstract 

> , 

In this paper, we give a complete analysis of an SIS epidemiological model in 
a population of varying size with two dissimilar groups of infective individuals. It 
is mainly based on the discussion of the existence and stability of equilibria of the 
proportions system and the result is in terms of a threshold parameter which governs 
j the stability of the disease free equilibrium. 
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1 Introduction 

The social mixing structure of a population or a group of interacting populations play a 
crucial role in the dynamics of a disease transmission. (See |§ and references therein.) In 
almost all attempts to combine epidemiological data with mathematical modeling, there 
has been a recognition of the need to consider the structure of social interactions among 
the individuals in the populations. (See [1 1] and references therein.) Many authors have 
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considered the multigroup models in which heterogenus subpopulations may participate 
to the epidemic process with different parameters 0. For SIS type models, a rather 
complete analysis of existence and global stability of a nontrivial epidemic state has been 
carried out by Lajmanovich and York JF]. In their work, the size of each subpopulation 
is assumed to be constant. 

A famous example for these subpopulations is the core group, i.e. the highly sexually 
active subgroups |J. It has become increasingly clear that the transmission within and 
among core subgroups is an important factor in the transmission of HIV/ AIDS ||. In 
order to consider the core group in an SIS epidemiological model, we divide the population 
into two subgroups each of them consists of susceptible and infective individuals. One of 
these subpopulations can be viewed as the core group. In this paper we consider another 
type of core group that is post-infection core group , i.e. individuals that become part 
of the core group after being infected. This hypothesis is plausible for a contagious and 
fatal disease like AIDS. From the psychological perspective, this group may be classified 
as violent. 

In this paper we examine an SIS model of disease transmission in a population of 
varying size with two dissimilar groups of infective individuals. One of these groups 
can be viewed as the post-infection core group . We also assume that the birth rate of 
susceptibles may be more than that of infectives. This is similar to the demographic 
assumption in M. This paper is mainly based on the discussion of the existence and 
stability of equilibria of the proportions system. First of all, in the next section, we 
introduce the model and some concepts of ODE's related to the system. In Section 3, we 
present some basic results concerning the nonexistence of certain types of solutions. In 
Section 4, we give a complete global analysis of the proportions system which is reduced 
to a planar system. The result is in terms of a threshold parameter which governs the 
stability of the disease free equilibrium. 



2 



2 The Model 



In order to derive our model, we divide the population into three groups: Susceptibles, 
S, and two groups of infectives, Ii and I 2 . We set N = S + Ii + I 2 which is the total size 
of the population and we use the following parameters which are assumed to be positive 
unless otherwise specified: 

b: per capita birth rate of susceptibles, 

b\\ per capita birth rate of infectives which is assumed to be < b, 

d: per capita disease free death rate, 

e: excess per capita death rate of infectives, 

Ai: effective per capita contact rate of I\, 

\ 2 : effective per capita contact rate of I 2 , 

71 : per capita recovery rate of I\, 

72 : per capita recovery rate of I 2 . 
We also assume that the susceptible individuals which have been infected, enter to the 
group I\ and I 2 of proportions p and q respectively, hence p + q = 1. 

The above hypotheses leads to the following system of differential equations in R^, 
where ""' denotes the derivatives with respect to t, the time, 



where b 2 = b — b\ and -jr- is of the proportionate (or random) mixing type , Q . By 
adding the above three equations, the total population equation is 




blN + (b 2 - d)S + 7i/i + l2 I 2 -Xi'-f- \ 2 '-f, 



(2-1) 
(2-2) 
(2-3) 



N' = (6i - d)N + b 2 S - e(h + I 2 ) 
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Setting s = jf,h = jf an< ^ ^2 = 77, we arrive at the following system of equations: 

s' = 6i(l - s) + b 2 s(l -s)+ 7l ii + l2 i 2 + (e - Ai)iis + (e - A 2 )i 2 s, (2 - 1)' 
= ps(Ai^i + A 2 i 2 ) + + - (h + e + ^y 1 )i 1 - b 2 si 1 , (2 - 2)' 

i' 2 = qs(X 1 i 1 + A 2 i 2 ) + ^2(^1 + — (h + e + ^ 2 )i 2 - & 2 si 2 - (2 - 3)' 

In order to determine the asymptotic behaviour of the solutions of this system of 
equations, we need the following concepts of ODE 's related to our system. 

Given an autonomous system of ordinary differential equations in IR n , 

§ = (2-5) 

we will denote by x.t the value of the solution of this system at time t, that is x initially. 
For V C W n , JCI, we let V.J = {x.t : x E V,t E J}. The set V is called positively 
invariant if VM + = V. For 7 CR" the cu-limit (resp. the a-limit) set of Y is defined to 
be the maximal invariant set in the closure of Y.[0, 00) (resp. Y.(— 00, 0]). We say that 
j(t) is an orbit running from x to x\ if lim j(t) = Xo and lim j(t) = x\. These two 
points must be equilibria and such an orbit is called heteroclinic orbit. When x\ coincides 
with x , it is called a homoclinic orbit. A closed curve connecting several equilibria whose 
segments between successive equilibria are heteroclinic orbits is called a phase polygon. 
By a sink we mean an equilibrium at which all the eigenvalues of the linearized system have 
negative real parts. Such a point is called a source if all of these eigenvalues have positive 
real parts. If some of these eigenvalues have positive real parts and the others negative 
real parts, then the equilibrium is called a saddle point and it is called nondegenerate if 
all of these eigenvalues are nonzero. 



3 Some Basic Results 

We start our analysis with some basic results about the system (2 — 1)' — (2 — 3)'. If we 
set £ = s + i\ + i 2 , then E' = (1 — E)(&i + b 2 s — ei\ — ei 2 ). Therefore the plane = 1 is 
invariant. We consider the feasibility region 

D = {(s, i 1 , i 2 ) : s + i x + %i — 1, s > 0, i 1 > 0, i 2 > 0} 
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which is a triangle and on its sides we have: 

s = s' = bi + 7i«'i + 7222, 
i\ = =^> z' x = pA 2 s2 2 , 
«2 = =>- z 2 = gAisii. 

It follows that .D is positively invariant and the disease free equilibrium (1,0,0) is the 
only rest point on dD, the boundary of D. Indeed our vector field points inward on 
dD — {(1,0,0)}. So every solution of the system (2 — 1)' — (2 — 3)' which starts in 

o 

dD — {(1, 0, 0)}, immediately gets into D, the interior of D. 

From now on, we examine the dynamics of this system in the feasibility region D. The 
following theorem is a modification of Theorem 4.1 in concerning the nonexistence of 
certain types of solutions. 



Theorem 3.1. Let / be a smooth vector field in R 3 and j(t) be a closed piecewise 
smooth curve which is the boundary of an orientable smooth surface S C M. 3 . Suppose 
g : U — > M 3 is defined and is smooth in a neighborhood U of S. Moreover it satisfies 
9{l(f))-f{l{t)) — an d (curl g) -n < 0, where n is the unit normal to S. Then 7 is not a 
finite union of the orbits of the system (2-5). 



In order to apply the above theorem, we define g = g\ + g 2 + #3 by 



0. 



Hi2 



Hi2 



n /i(a,»a) 
f2{s,h) /i(s,*i) 







where /1, /2 and ^3 deduced by S = 1 on the right hand side of (2 — 1)', (2 — 2)' and 
(2 — 3)' respectively. Now after some computations we get 



(curl g). (1,1,1) 



p\ 2 : q\ 1 : bi + 71 61 + 72 



llS* 



Corollary 3.2. The system (2 — 1)' — (2 — 3)' has no periodic orbits, homoclinic orbits 

o 

or phase polygons in D- 
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Proof. We use Theorem 3.1. for / = (fa, fa, fa). Here we have g.f = and (curl <?).(1, 1, 1) < 
in D. □ 



Lemma 3.3. The cj-limit set of each orbit of the system (2 — 1)' — (2 — 3)' with initial 
point in D is a rest point. 

o 

Proof. Suppose the contrary, then the cu-limit set has a regular point in D- Let x be 
such a point and h be its first return map. For a point y near x on the transversal, let V 
be the region surrounded by the orbit 7 from y to h(y) and the segment between them. 
This region is known as Bendixon sack. (See Figure 3.1.) 

Now by Stokes' theorem 

/ / (curl g). (1,1, l)da = [g.fdt+ [ g(ty + (1 - t)h(y)).(y - h(y))dt. 

J JV J-y JO 

Since g.f = and h(x) = x, the right hand side of the above equality tends to zero when 
y tends to x. But the left hand side tends to the integral over the region bounded by the 

o 

tu-limit set. This is a contradiction since (curl g).(l, 1, 1) < in D- □ 




Fig. 3.1. The Bendixon sack. 



Remark 3.4. When the u;-limit set lies in D the above result is easily concluded by the 



generalized Poincare-Bendixon theorem [K| and Corollary 4.2. Similarly if the a-limit 
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o 

set of an orbit of the system (2 — 1)' — (2 — 3)' lies in D, it must be a single point. 



4 The Planar System 

Using the equality s + i± + i 2 — 1, we see that our system is essentially two dimensional. 
Thus we can eliminate one of the variables, say s, to arrive at the following quadratic 
planar system 

i'i = (pAi - b-e - 71)21 + pX 2 i 2 + (h + 12) ((h + e - pX^h -pX 2 i 2 ), (3 - 1) 
i' 2 = q\ l i l + (q\ 2 - b - e - ^ 2 )i 2 + (h + ^)((&2 + e - gA 2 )i 2 - qX x ii). (3 - 2) 

The dynamics of the system (2 — 1)' — (2 — 3)' on D is equivalent to the dynamics of this 
planar system in the positively invariant region 

Di = {(M 2 )|ii > 0,i 2 > 0,h+i 2 < 1}. 

The matrix of the linearization of the system (3 — 1), (3 — 2) at the origin is: 

p\ 1 - b - e - 71 pX 2 

qXi qX 2 - b- e — 72 

with det C = (6 + e + 7i)(6 + e + 72) - pAi(6 + £ + 72) - qX x {b + e + 71). We set 
Rq = + Hence if i? < 1, then det C > and trace C < and if R > 1 

then det C < 0. Thus we have proved the following lemma. 

Lemma 4.1. Let R be the above threshold. Then the origin is a sink (resp. a saddle) 
for the system (3 — 1), (3 — 2) whenever R < 1 (resp. R > 1). 

Lemma 4.2. The trace of the linearization of the system (3 — 1), (3 — 2) at a rest point 

o 

in Di is negative. 
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Proof. We compute the trace at a rest point in D\. 

|^ = p\i - b - e - 71 + (b 2 + e - p\i - p\ 2 )i 2 + 2(6 2 + £ - pXi)h, 
§| = pA 2 - 6 - e - 72 + (6 2 + e - q\ 2 - gAi)ii + 2(6 2 + e - gA 2 )i 2 . 

From i[ = and i' 2 = 0, we get 

§| = -pA 2 | + P A 2 | + (b 2 + e- p\ 1 )i 1 = -pA 2 f(l - i 2 ) + (6a + £ - pAi)ii, 
g = -gAxg + gAx| + (6 2 + e - q\ 2 )i 2 = -gA^(l - ix) + (fo 2 + £ - gA 2 )* 2 . 

Using the equality s + i± + i 2 — 1, we obtain 

<9« 2 /, x \ /, x \ x ^ 2 s , ixs 
77- + 77- = (»2 + e - Ax)«x + {b 2 + e- \ 2 )i 2 - pA 2 - qXi — . 

OTx C/Z 2 lx 2 2 

Now from (2 — 1)' we have 

s' = 6x(«i + ^2) + &2s(^i + ^2) + 7i«i + 72^2 + (e - Ax)ix« + (e — A 2 )? 2 s = 0. 
Thus (6 2 + £ - Ax)ixs + (6 2 + £ - A 2 )i 2 s < and it follows that §| + §| < 0. □ 

The following two corollaries are immediate results of the above lemma. 

o 

Corollary 4.3. The system (3 — 1), (3 — 2) has no source in Di- 

o 

Corollary 4.4. Every nondegenerate rest point of the system (3 — 1), (3 — 2) in D± is 
hyperbolic. 

Remark 4.5. A nondegenerate rest point of the system (3 — 1), (3 — 2) is obtained by 
a transversal intersection of the two conic sections i[ = and i' 2 = 0. 

o 

Proposition 4.6. There is at most one rest point in Di for the system (3 — 1), (3 — 2). 
Moreover such a rest point is always hyperbolic. 
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Proof: From the equilibrium conditions i[ — i' 2 — 0, we get the following equation which 
is homogeneous with respect to i\ and i 2 of second order. 

((pAi - 6 - e - 7i)ii + pA 2 i 2 )((6 2 + £ - <?A 2 )2 2 - gAiii) = 

(gAiii + (gA 2 -b-e - 7 2 )i 2 )((6 2 + £ - pAi)ii - pA 2 i 2 ). 
This equality can be written as 

gAi(6i + 71)1? + (*)iii 2 - pA 2 (&i + 72)^ = 

where (*) is a statement in terms of the involved parameters. The set of all roots of 
this quadratic equation consists of two lines through the origin in the (ii,i 2 ) plane. One 
of these lines has negative slope and meets D 1 only at the origin. Thus the other line 
contains all rest points of the system (3 — 1), (3 — 2) in D\. Since each line contains at most 
two rest points of a quadratic planar system and this line contains the origin, it follows 

o 

that Di contains at most one rest point. This rest point is obtained by a transversal 
intersection of this line and each of the conic sections i[ = or i' 2 = 0. It is easy to 
see that at this rest point, these two conic sections must intersect transversally. Now by 
Remark 4.5, this rest point is nondegenerate and by Corollary 4.4, it must be hyperbolic. 
□ 

Remark 4.7. We have indeed shown that all rest points of the system (3 — 1), (3 — 2) 
which are not more than three points, are nondegenerate, except the origin in the case 
#0 = 1. 

Now we are ready to prove our main result about the dynamics of the system (2 — 
1)'- (2-3)' in D. 

Theorem 4.8. (i) If R < 1, then (1, 0, 0) is a global attractor in D 

o 

(ii) If i? > 1, then there exists a unique rest point (an endemic equilibrium) in D which 
attracts D - {(1,0,0)}. 

Proof. When R < 1, the origin is a sink for the planar system (3 — 1), (3 — 2). If there 

o 

exists another rest point in Di for this system, it must be unique and hyperbolic. By 
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Corollary 4.3 it cannot be a source. If it is a sink, then we will have two sinks in D\. 
The basins of attraction of these two points are open and by Lemma 3.3, D\ is the union 
of these two open subsets. This contradicts the connectedness of D\ and shows that it 

o 

cannot be a sink. Now suppose that there is a saddle point in D\- By Corollary 3.2, 

o 

there is no homoclinic orbit in D\- Hence the origin attracts the unstable manifold of the 
saddle point. Now the region bounded by the unstable manifold contains some part of 
the stable manifold of of the saddle point. Thus the a-limit set of this part of the stable 
manifold is a rest point in Di by Remark 3.4. This is a contradiction since this rest point 
can be neither the origin nor the saddle point. Therefore the origin is the only rest point 
in D\ for the system (3 — 1), (3 — 2) and by Lemma 2.3, it is the u;-limit set of all points 

of £>!. 

The above fact is still valid for the limiting case, R = 1. To see this, suppose that 

o 

there exists another rest point in Di, then it must be hyperbolic and belong to D\- Thus 

o 

it remains in Di when the involved parameters are slightly changed to get Ro < 1 which 
contradicts the above result. This finishes the proof of (i). 

Now suppose R > 1. Then the origin is a saddle point for the planar system (3 — 

o 

1), (3 — 2). Thus by Lemma 2.3, there must be some rest point in Di- Since such a rest 
point is unique and hyperbolic, it must be a sink and attract all points of D\ except the 
stable manifold of the origin. We claim that the stable manifold meets D 1 only at the 
origin. To see this notice that some part of the unstable manifold of the origin must be 
outside of D\ for its right angle. Since D\ is positively invariant, the stable manifold 

o 

does not intersect Di- Moreover, the vector field points inward on dD\ — {(0,0)}. Thus 
the stable manifold of the origin does not intersect dD\ — {(0,0)} either. This shows 
that the origin cannot attract any point of D\ — {(0,0)}. Thus the unique sink attracts 

o o 

Di —{(0, 0)}. It means that there is a unique rest point in D for the system (2—1)'— (2— 3)' 
which attracts D — {(1, 0, 0)}. □ 

Remark 4.9. In the above argument, in order to prove the global asymptotic stability 

o 

of the endemic equilibrium (i.e. the unique rest point in Di), we showed that the stable 

o 

manifold of the origin cannot intersect Di- It is a special case of the following fact. Let 
X be a smooth vector field on a smooth manifold M and D C M is a positively invariant 
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region with a saddle point on dD. If the unstable manifold of this saddle point contains a 

o 

point of [M — D)°, then its stable manifold cannot intersect D- In order to prove it, one 
can follow our proof in the special case and observe that this is a direct consequence of 
the Hartman-Grobman theorem. However, it is obvious by the Inclination Lemma [Ibl . 



Remark 4.10. If we consider I2 as the post-infection core group, then the threshold 
-^o — b+e+j! + b+e+72 c l ear ly shows the effect of this group on the epidemics process. 
Although the probability q is a small number, the fact A 2 >> Ai causes the term q\ 2 to 
be significant. 
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